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Abstract 

We investigate the classical ground state configurations of a collection of skyrmions in the limit 
of vanishing Lande g factor. We show that in this limit, the skyrmions which are large and 
overlapping, behave qualitatively differently from small skyrmions. We investigate the system in 
various regimes of spin stiffness. We find that in all regimes, at = 0, the ground state configuration 
is a rectangular skyrmion lattice with Neel orientational ordering. The charge distribution has a 
higher symmetry of a face-centered rectangular lattice and the z-component of the spin polarization 
is zero. We also present general argument for the vanishing of spin polarization at g( = at finite 
temperature and in the presence of disorder and discuss some experimental consequences of our 
results. 
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I. INTRODUCTION 



The lowest energy charged excitations about the z/ = 1 ferromagnetic quantum Hall 
ground state are skyrmions. They are topological objects in which the spin gradually twists 
over an extended region. Their spin is greater than 1/2 and they carry an electric charge 
of ±e. Skyrmions have a topological charge which equals their electric charge at u = li^. 
Skyrmionic excitations are favoured over single particle excitations when the Lande g-factor 
is small i.e g — >■ 0^*^. Skyrmions which are produced around the u = 1 ground state have 
been experimentally seen by in OPNMR and optical magneto-absorption experiments on 
Ga nuclei in an electron doped multiple quantum well structure^. A sharp fall in the spin 
polarisation of the 2D electron gas is seen on either side of = 1, which indicates that a 
number of spins are being flipped by the addition or removal of a single electron, in contrast 
to single particle excitations where no extra spin is flipped. Small skyrmions have also 
been observed as low energy excitations about the u = 1 state in tilted field measurements 
reported in^. The above experiments observe skyrmions where the number of flipped spins 
varies between 3 and 7. Low g experiments observe large skyrmions with a large number of 
flipped spins (50 and 36 respectively)^'^. In these experiments a low value of g is obtained 
either by suitably altering the stoichiometry while growing the sample or by application of 
hydrostatic pressure. 

The low temperature ground state of a system of interacting skyrmions is expected to be 
a crystalline lattice^SJiiSii^. Earlier calculations performed in the mean field limit by Brey 
et. alM suggest that the ground state of a system of a skyrmions is a square Skyrme crystal 
with an anti-ferromagnetic order in the planar component of the spin of the skyrmion. Later 
studies by Green et. al. , Rao et. al. and Timm et. al. analyze the classical ground state of 
the skyrmion lattice as — 1 using the non-linear sigma model (NLSM) ^^i^^i^^ . Green et. al 
assume that the ground state is a triangular lattice with spiral orderingii-. Rao et. al. report 
triangle square —>■ triangle transitions in the skyrmion lattice as a function of the filling 
factor near u = 1 and at T=C^^. Timm et. al. study a system of well separated skyrmions as 
described by an anti-ferromagnetic XY model and propose a T=0 phase diagram in which 
the triangle and square phases are separated by Neel ordered, face-centered rectangular 
phases^^. The classical and quantum phase transitions occurring in the Skyrme crystal has 
been studied using Hartree-Fock calculations by Cote et. a/.^^. Moon et. al. study the 
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renormalization effect of skyrmions and the thermodynamics of a skyrmion system, without 
assuming apriori lattice structures'^. 

At non-zero g, the charge and spin densities of skyrmions are peaked about its north 
pole if the vacuum spin polarization direction is taken to be the south pole. They fall 
off exponentially from their central value. The corresponding size is determined by the 
competition between the Coulomb repulsive energy and the Zeeman energy. The former 
favours an infinite size (power law fall-off) and the latter a zero size. The XY component of 
the skyrmion spin density assumes a "hedgehog" configuration. These configurations can be 
characterized by the orientation angle subtended between the position vector (with respect 
to the center of the skyrmion) and the XY component of the spin density at that point. 
The finite sized skyrmion therefore behaves like an oriented, charged disc. It is energetically 
favourable for nearby skyrmions to have opposite orientations. These considerations lead to 
small skyrmions being modeled by discs carrying an XY spin with coulomb repulsion and 
short range anti-ferromagnetic interactionsti^. The coulomb repulsion leads to a triangular 
lattice in the dilute limit and the XY antiferromagnetism in the orientational degree of 
freedom to the square lattice as the inter-skyrmion spacing decreases. 

In this paper we show that this picture breaks down when the skyrmionic size is larger 
than the inter-skyrmion spacing. For instance, when g = there is no exponential damping 
and we show that the charge density is symmetric in the neighbourhood of the north and 
south poles. The oriented disc picture is then clearly invalid. We study the ground state 
configuration of a collection of skyrmions at g = and g = 0.05go, {go is the physical 
value of the Lande g factor). The energy of the system at (7 ~ is described well by 
the Non-linear Sigma Model energy functionatf^. We minimize this energy functional in 
three different regimes: (1) High stiffness regime in which the gradient term dominates the 
energy functional, (2) Low stiffness regime in which the Coulomb term dominates and (3) 
Intermediate stiffness in which both terms in the energy functional compete to generate a 
suitable ground state. Finally, we discuss the behaviour of the spin polarisation and its 
experimental consequences. 

The rest of the paper is organised as follows. Section II, describes the NLSM energy 
functional. Section III discusses the ansatz used to solve the energy functional and the 
symmetries of that ansatz. In Section IV we discuss the minimum energy configuration in 
the high stiffness ,low stiffness and intermediate stiffness regimes. Section V we describe 
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a general argument for the vanishing of the spin polarization at g = and discuss the 
experimental consequences. Finally, the conclusions of this paper and a summary of the 
results is presented in Section VI. 



II. MODEL AND NLSM ENERGY FUNCTIONAL 

We consider skyrmionic excitations about the u = 1 ground state. The number density 
of skyrmions is given by Usky = ^v^^c, where ric is the carrier density. In our calculations we 
use the carrier density, Uc = 1.5 x 10^^ cm~^ and change u by tilting the magnetic field. We 
localize the skyrmion centers on the lattice points of a Bravais lattice shown in Fig. 1. We 
choose a skyrmion density (usky) of 1 skyrmion per unit cell (i.e one skyrmion per lattice 
point). 

We measure all lengths in units of the magnetic length l^, = and all energies in terms 
of the cyclotron energy huOc-, where uo^ = {eB/m*c), m* is the effective mass of the electron. 
The area of the unit cell is fixed by the filling factor and is given by 

= ^ = ^ (2.1) 

The local spin polarisation which is represented by unit vector field n(x) is stereographi- 
cally projected onto the complex plane by the transformation w = cot(6'/2)e*'^, where 6 and 
are the polar angles of the spin vector n(x). In the rest of the paper we will work with 
the planar spin variable w. 

The topological charge density is given by 

p(x) = -^eijn.(9in x djii) 

e.jdiwdjw ^2 2) 



2Tci{l + wwy 

The topological charge, (^(x) = /□ d'^x p(x) = 1 (/□ denotes integration over the unit cell). 

The low energy, long wavelength excitations about the u = 1 ground state are described 
by the NLSM. The NLSM energy functional has to be minimized for different filling factors 
to get the minimum energy configurations of this lattice. The NLSM energy functional with 
Zeeman and Coulomb interactions 13^*^^*^^: 

E = Egrn + -Ez + Ecoul (2.3) 
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The gradient or the spin exchange term proportional to / d?x\ c}jn(x) = x, y), is 
calculated as 



E,rn =i L rf2^ V-----";^2r'-^ (2.4) 



7 f {dr^wd^w + dywd yw) 
2 In ^ {l + wwY 
where 7 = The gradient energy density is Egm/An , where Au is the area of the 

unit cell. This term alone is the pure NLSM and it has exact scale invariant solutionsi^. 

The Zeeman term is proportional to the z-component of the total spin i.e. to 
^ /□ (Px ^^^^ (where is the z-component of n(x)). In our units the average number 
of electrons is The z-component of the total spin is 

U f 1V1V 

(Total spin)z = — / d'^x-. _ , (2.5) 

2n Jn [1 + WW) 

Therefore, 

E.^s'ffj^x-f^^ (2.6) 

271 Jn [1 + WW) 

where g* = ^f^- The Zeeman energy density is therefore E^/ Aq, where An is the area of 
the unit cell. 

The Coulomb energy density term is a term of the form 

Ecoui='^^[ p(x)^^p(y) (2.7) 
2 Atot J^,y I X - y I 

where e* = {e'^/Klc){l/fiuJc) and Atot is the total area of the lattice. 

The Coulomb term arises because the electric charge density is proportional to the topo- 
logical charge density. Since the topological charge density explicitly appears in the above 
expression the spin orientation gets automatically tied to the Coulomb energy. The four 
dimensional integral in the Coulomb term can be converted to a sum over the reciprocal 
lattice. 

If 1 

Ecoui = / 

An "'x.yen | X — y + K | 

After doing the summation over {R} and converting the integral over g to a sum we get 

Ecoui = ^ E I KGr) (2.9) 

where p{Gr) = /^g^ p{'x)e^^^-^, Gr lies in the reciprocal lattice and An is the area of the 
unit cell. 
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III. ANSATZ AND SYMMETRIES 



In this section we will describe and motivate the ansatz we will use to minimise the 
energy functional discussed above. Let us first consider the model without the Coulomb and 
Zeeman interactions. In this case it is known^'^ that any meromorphic function w{z) is an 
exact solution of the equations of motion. Since the positions of the poles and zeros define a 
meromorphic function upto a multiplicative constant, the most general skyrmion solution 
can be written as, 



Wo{z) = Woo /. — (3-lJ 



nli(^-R„) 

Woo specifies the direction of the spin polarization at \z\ oo. Equation (jH.lj) can be 
rewritten as, 

n=l \^ ^n) 

\n are complex numbers which are given by the solution of the N linear equations, 

^oo + Ep-^ = (3.3) 

If we choose Woo = 0, corresponding to the spin polarization at infinity pointing in the 
—z direction, and take the N ^ oo limit, equation ()3.1|) is written as, 

-oW^E^ (3.4) 

This motivates an ansatz of the form 

w{z) = Wo{z)F{z) (3.5) 

for the fully interacting model. The skyrmion lattice is obtained by letting {R} go over 
the position vectors of a lattice. For a Bravais lattice with basis vectors el and e2, R = 
nel + me2, where n and m are integers. Further when |el| = |e2| = e these basis vectors 
describe a face-centered rectangular lattice as shown in Fig. 1. They subtend an angle 6 at 
the origin. When 6 = 271/3 the lattice is a triangular lattice. When 9 = 71/2 the lattice is 
rectangular with a and h as the length of its sides, a/6 is the aspect ratio of the rectangle. 
a/h = 1 corresponds to the special case of a square lattice. F{z) is a smooth, non-analytic 
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function with the periodicity of the lattice. For the rectangular lattice we take it to be of 
the form, 

F{z) = 

g ^sp'^sp 

Tnp = \J sm^iiix / a) + sin^(7r?//6) 
Tsp = \J cos^^-Kx/ a) + cos^(7ry/6) 

The spin configuration corresponding to this ansatz is quasi-periodic. We have, 

w{z + K) = e'^^-^wi^z) (3.6) 

i.e. the spin rotates about the 2;— axis by and angle q.R when translated by R. However 
the energy has the periodicity of the lattice (since all the energies depend on ww. As shown 
in earlier work^^, w{z) describes a configuration with skyrmion number one per lattice site. 
The mapping w = cot{6 /2)e^'^ implies that the anti-podal points = and 9 = n correspond 
to the poles and zeroes of wo{z). In the rest of the paper, we will refer to the poles as the 
"north poles" and the zeroes as the "south poles" respectively. 

If we take Kgp = and 1/nnp << ci,b, then the ansatz in equation ()3.5|) is essentially the 
same as the one used in^^. The charge and spin densities are then concentrated about the 
north poles. In this regime, each skyrmion can be visualised as a charged, oriented disc. 
The orientation corresponding to the angle the XY component of the spin density at any 
point makes with the vector joining that point with the center of the disc. The ansatz in 
equation ()3.5p then corresponds to a collection of well separated skyrmions centered at the 
lattice sites (see Fig. 1). The vector q can be interpreted as specifying a spiral orientation 
order for these skyrmions. The orientation rotates by q.R when we move by R. This is the 
picture used in reference^^ and the classical phase diagram can be qualitatively understood 
in terms of the long range coulomb repulsion between the discs (which favours a triangular 
lattice) and a short range antiferro-orientational interactions (which prefer an unfrustrated 
square lattice). When g = —0.44, the value for GaAs at low pressures, previous theoretical 
and experimental work show that the skyrmions have a spin s ~ 5. The spin can be related 
to the rms charge radius as (r^) = 2s (see Appendix 2). Skyrmions can be called small when 
this radius is smaller than half the lattice spacing ^J^^- Skyrmions of spin ^ 5 have a size 
of ~ 3. At z/ = 0.9, the lattice spacing is around 8. The skyrmions are thus well separated at 
V > 0.9. The system of skyrmions in this regime is well described by the oriented, charged 
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disc picture. As g is decreased, the skyrmions become larger and, as we will show below, 
this picture breaks down. 

We show plots of the spin density and charge density in a lattice of well separated 
skyrmions. Figures 2. and 3. show the XY spin configuration and the charge density 
in a triangular lattice of well separated skyrmions^A. The spins are ordered into a triangular 
lattice (with ABC sub-lattice structure) with q = (271/3, 27r/3). Figures 4. and 5. show the 
XY spin configuration and the charge density distribution for a square lattice of well sepa- 
rated skyrmions. The spins have Neel order (with AB sub-lattice structure) with q = (tt, vr). 
We notice that the charge is peaked only around the "north poles" (the lattice points). 

Let us now consider the extreme case i.e. = 0. We first look at the high stiffness limit 
i.e. when the coefficient of the gradient term, 7 is large. The gradient energy then dominates 
the NLSM energy functional. This energy is exactly minimised by any analytic ansatz. So 
in this regime, we will be putting = = 0. The ansatz is then w{z) = Wo{z) with q 
and A as the variational parameters. The energy minimization procedure and results will 
be described in the next section. Here we discuss some of the important features of the spin 
and charge density distributions corresponding to this ansatz. 

Since there is exactly one skyrmion per unit cell, there is one north pole and one south 
pole per unit cell. From the form of wo{z) given in equation ()3.4|) . it is obvious that the 
north poles are at {R}. As we will now show, for a given lattice, q determines the positions 
of the south poles. First consider a system of two skyrmions located at R and — R. This 
configuration is described by, 

„iq.R „— jq.R 

M^) = A(^ + ^) (3.7) 
z — K z + lt 

The zero of this ansatz is located at — 'iRtan(q.R). Thus given the positions of the two 
skyrmions, the position of the zero is determined by q. 

We will now illustrate that the position of the zero is determined by the value of q by 
considering the case of a face-centered rectangular lattice. These lattices are described by 
the basis vectors el and e2. The unit cell is also shown. The center of the unit cell is located 
at Zq = (el + e2)/2. Let us first consider Neel order q = (vr, tt). Now, 

gi7r(n+m) 

= ^- (n-l/2)el + (m-l/2)e2 ^^'^^ 
Reflecting every lattice point about the point Zq by the transformation n = — n' — 1 and 



m = m' — 1, we can show that the above equation becomes 

gi7r(n'+m') 



{-n' + l/2)el + (-m' + l/2)e2 
= -w{zo) (3.9) 

which imphes w{zq) — 0. Hence the point zq which is the center of the unit cell of the 
Bravais lattice is the location of the zero of the ansatz or the south pole. We note here that 
the point zq — (el + e2)/2 is the south pole independent of specific values of el and e2 as 
long as there is Neel order. Hence in a triangular lattice and a rectangular lattice the south 
pole will lie at the center of the unit cell when these lattices are Neel ordered. 

Now consider a triangular lattice with q = (27r/3, 27r/3) ordering. We will now show that 
the position of the zero lies at the point cq — (el — e2)/3. A rotation of a lattice point 
R about this point through an angle of 27r/3 takes it to a point R' which also lies on the 
lattice. We know, 

gi27r/3(?i+m) 

w{co) = Snm (g^_R) (^-lO) 

On multiplying e*^"^/^ to the numerator and denominator of w{co) it becomes 

gi27r/3(n'+m') 

w{co) = e'^^/'j:n'm'-, ^ (3.11) 

(Co - R'j 

Since R' is also a Bravais lattice vector 

wico)=e''^/'w{co). (3.12) 

Hence, we find that cq = (el — e2)/3 is the south pole of the ansatz for a triangular lattice 
with (27r/3, 27r/3) order. 

When the ordering on a triangular lattice changes from Neel (tt, tt) to (27r/3, 27r/3), the 
position of the south pole moves from zq — (el + e2) /2 to the point Cq = (el — e2) /3. Hence 
we see that the position of the south pole of the ansatz on a fixed lattice is controlled by 
the parameter q. In the rest of this section we will consider a rectangular lattice with Neel 
order and an aspect ratio of a/6. The south poles of the ansatz are located at the centers of 
the unit cells. 

We will now study the properties of the ansatz on a rectangular lattice with Neel order. 
We show that there exists a special value of A = at which the charge and spin densities 
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are symmetric about the north poles and south poles. This is true if the ansatz obeys the 
condition 

w(z+ (el + e2)/2) = at A = A, (3.13) 

W{Z) 

and hence the charge density and spin polarisation density are 

pH = ^TiW = /^(V^) (3.14) 



sjw) = f = -sJl/w) 

Thus, if Eq. 13.131 is true, the charge density is symmetric near the north and south poles 
and the spin is also symmetric but in opposite directions near these poles. 

We impose the condition in Eq. ()3.13|1 close to the north and south poles and obtain an 
equation for A^ as follows. Expanding w{z = Znp + e) about the north pole (znp) we find that 
to leading order w{z) = ^. Expanding about the position of the south pole zq = (el + e2)/2 
we get 

w{z) = w{zo) + edzw{z) 

A(-l)'^l+'^2 

= '^"^'"^(^o-(r.iel + n2e2))2- ^^'^^^ 
The condition in Eq. 13.131 means that As obeys 

Af = ^"^''^nzo-(niel + n2e2))2- ^^-^^^ 

This was solved for numerically for As and the condition in Eq. 13.131 was checked for all 
z. For example, for a Neel ordered rectangular lattice with an aspect ratio of y^j As was 
numerically calculated to be equal to 8.42. Even though the As was obtained by imposing 
Eq. 13.131 at points close to the north and south poles, the numerical calculation described 
above shows that Eq. 13.131 is satisfied at all z. Thus at Ag the charge and spin densities are 
symmetric about the north and south poles. 

Having shown that at a particular value of A = As the charge and spin densities are 
symmetric around the north and south poles of a rectangular lattice with Neel order, we 
motivate some more properties of this configuration which will be numerically verified in 
the next section. For such a rectangular lattice with an aspect ratio of v^, the north poles 
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and south poles form a smaller triangular lattice. At the value A = A^, the charge and spin 
density also have the symmetry of this smaller triangular lattice since they are symmetric 
about the north and south poles (see fig. 3). Since such a configuration corresponds to a 
charge distribution with the triangular lattice symmetry we expect that the Coulomb energy 
will be minimised for such a distribution. For the special case of the square lattice with aspect 
ratio 1, the charge and spin densities also have the symmetry of a smaller square lattice at 
the value A^. This is our motivation for minimising the variational ansatz for w{z) with the 
variational parameters as 6, q, A and aspect ratio. We do this calculation numerically to 
find the optimum value for these parameters in the lowest energy configuration. 

In order to contrast the behaviour of spin and charge density for overlapping, large 
skyrmions with that of well separated, small skyrmions we plot the XY spin configurations 
and charge density profiles for a lattice of overlapping skyrmions. Figures 6. and 7. s how 
these plots for a rectangular lattice with an aspect ratio of ^/S. We see that the XY spin 
configuration is very different from that shown in Figure 2. Also, the charge orders into a 
honeycomb lattice which has the symmetry of a triangular lattice. Figures. 8. and 9. show 
the plots for a square lattice of overlapping skyrmions. 

We note here that the above features of the analytic ansatz remain true for a general 
ansatz w{z) = wo{z)F{z) as long as F{z) is a function that respects the north pole-south 
pole symmetry. We discuss below a certain type of such functions and also conditions under 
which the north pole- south pole symmetry is broken. 

In the low and intermediate stiffness regimes the gradient term does not dominate the 
energy functional anymore. Hence the requirement of the analyticity of the ansatz is not 
valid any longer. The ansatz can have an exponential cut-off via the function F{z). In view 
of the north pole-south pole symmetry (which is still valid at (7 = 0), we choose to damp 
the ansatz around the north and south pole in a symmetric way. Accordingly we choose, 
i^np = i^sp in Eq- ()3.6|1 . This form for F{z) preserves the north pole-south pole symmetry. 

Symmetry about the north and south poles can be broken in two ways. One way is by 
choosing the scale A to be different from the symmetric value A^ thus making the charge and 
spin densities asymmetric about the north and south poles. The other way is to use different 
exponential cut-off scales around both poles {i.e Unp 7^ ^^sp)- This can be accomplished by 

^(^) = (3-17) 

p n,2 ' sp 
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Such an asymmetry can be induced by a large value of g. In the extreme case when the 
cut-offs are such that the charge is peaked entirely around the north pole, the overlap is 
very small and the oriented disc picture is valid. 

IV. ENERGY MINIMISATION : GROUND STATE CONFIGURATIONS 

In this section we show that the special value at which the north and south poles 
are symmetric is the one at which the energy is minimised. We also show below that in 
the minimum energy configuration the charge density has the symmetry of a triangular 
lattice in the high stiffness limit and that of a square lattice in the low stiffness limit. This 
corresponds to the skyrmion lattice being rectangular with an aspect ratio of \/3 in the 
former case and square in the later case. In intermediate stiffness regimes, the aspect ratio 
interpolates between \/3 and 1. At the minimum, the spin density is peaked equally in 
opposite directions at the north and south poles which means (s^) = in the unit cell. The 
vanishing of the ^-component of the spin the unit cell at A = is a universal feature of our 
calculations. 

A. High Stiffness Regime 

In this regime at (7 = 0, the NLSM energy functional has contributions from the gradient 
term and the Coulomb term. Since the coefficient multiplying the gradient term, 7 (the spin 
stiffness), is large the gradient energy dominates the energy functional. We minimize the 
energy functional with the analytic ansatz defined in the previous section. As mentioned 
earlier, this ansatz minimizes the gradient term exactly and hence is a good choice in the 
high stiffness limit. For large, overlapping skyrmions, we expect that the charge density has 
a higher symmetry than the underlying lattice. Based on the arguments in the previous 
section, we can see that this higher symmetry involves the symmetric distribution of charge 
and spin around the north and south poles. The charge density has the symmetry of a 
face-centered rectangular lattice if the underlying skyrmion lattice is rectangular. 

We numerically minimize the NLSM energy functional with the ansatz for w{z) on a 
rectangular lattice with the quantities q,A and aspect ratio a/h as the variational parameters. 
Neel orientation order [i. e. q = (tt, tt)) is chosen by the minimum energy configuration. 
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This ensures that the south pole hes at the center of the rectangular unit cell. An aspect 
ratio of a/h — -\/3 is chosen by the minimum energy configuration (see Fig. 10). At these 
values for q , a/6 and A = Xmin the charge density is symmetrically distributed around the 
north and south poles of the distribution as shown in Fig. 11. Hence this value of A = Knin 
is also the special value A^ discussed in the previous section at which there is complete 
symmetry between the north and south poles. The charge density has the symmetry of a 
triangular lattice at this point. 

The spin density peaks equally in opposite directions at the north pole and south pole. 
Hence the most symmetric distribution corresponds to a situation in which the total z- 
componcnt of the spin is zero. 

Because of the absence of the competing Zeeman term, there is no other length scale 
competing with the length scale of the Coulomb energy. The scaling relation is Ecouii^, q) — 
jEcouii^: q) (see Appendix 1). As the filling factor is varied, the length scale corresponding 
to the Coulomb energy changes. Thus, the length scale A should scale with the fiUing 
factor (lattice spacing) but the ground state configuration does not change. If the energy 
is minimised at a certain filling factor (corresponding to lattice spacing ei) by a value Ai, 
then the energy at another filling factor (corresponding to another lattice spacing 62) is 
minimised by a value A2 such that Ai/ei = A2/e2 is satisfied. This scaling of Xmin with 
the lattice parameter e is shown as u is changed (see Fig. 12). The ground state remains a 
rectangular lattice with aspect ratio of VS with Neel order for all filling factors at g — 0. 

To conclude this subsection we give some details regarding the actual numerical min- 
imization procedure used and some checks that were done. The numerical routine used 
was a downhill simplex routine in the three parameter space of q,A and a/b. This routine 
minimized the energy functional with respect to the above parameters and returned the 
minimised value of energy at the appropriate values of these parameters. The minimization 
was carried out for many fiUing factors close to 1. As explained in Section II, different filling 
factors correspond to varying the area of the lattice. The numerical code was first checked 
to give the correct topological charge of 1 per unit cell (with an error of less that 0.1%). 
The vanishing of the 2;-component of the spin was numerically seen when A = Xmin, lending 
support for our expectation of complete north pole-south pole symmetry at this point. 
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B. Low and Intermediate Stiffness Regimes 

We now consider the opposite limit (zero stiffness limit) at g = for a system of overlap- 
ping skyrmions. In this limit only the Coulomb energy contributes to the energy functional 
aX g = and hence minimum energy configurations of the skyrmion lattice are determined 
by the minimum of the Coulomb energy. The Coulomb energy is minimized by a suitable 
configuration in which the charge is uniformly distributed. The ansatz used for the min- 
imisation need not be analytic any longer and can have an exponential cut-off. Keeping in 
mind the north and south pole symmetry we chose to damp the ansatz around the north 
and south poles in a symmetric way. This is done by choosing F{z) to be a function with 
the periodicity of the north pole-south pole structure in the unit cell. The form for F{z) is 
given in Eq. ()3.6|) with Knp = i^sp = 

The energy is minimized with this ansatz on a Bravais lattice with respect to q, A, n and 
aspect ratio a/b. We find that the ground state configuration is a square lattice with Neel 
order and Sz = (see Fig. 13). The type of lattice does not change with filling factor and 
the scaling of A and Coulomb energy with lattice size is observed. The charge density has 
the symmetry of a square lattice. The spin density is also peaked in opposite directions at 
the north and south poles. 

We have also checked whether an ansatz which describes a charge density with a different 
behaviour at small r will change the features described above. For this purpose we chose an 
ansatz 

w{z) = woiz)F{z) 

Fiz) = ^"""^ 1^ + '^'-^l (4.1) 

ensuring the symmetry of the north and south poles. When the minimization was done for 
this ansatz we find that Neel order and = are still preferred with a square lattice as the 
ground state at all filling factors. At filling factors below u = 0.93, we find that k,2 << i^i 
and the charge density has the same behaviour as shown in Fig. 13. Between filling factors 
0.93 and 0.98 we find that ^2 > ni and the charge gets distributed to the lattice points as 
well as the center of the unit cell. 

In the intermediate regime the stiffness is chosen to be at its physical value 7. All 
three terms (gradient energy, Zeeman energy and Coulomb energy) contribute to the energy 
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functional. The ansatz used in Eq. ()3.5p is used to determine the ground state configuration 
at g = 0. We find that the minimum energy configuration is a rectangular lattice with an 
aspect ratio which is ansatz and filling factor dependent. However the spin configuration 
still has Neel order and = for all these configurations. 



V. EXPERIMENTAL CONSEQUENCES 

Our calculations predict two universal features of the g = system, (i) The skyrmion 
lattice is rectangular with Neel orientational order, (ii) the spin polarization is exactly zero. 
It is difficult to directly experimentally observe skyrmion lattices and no such experiments 
have been done so far. However the spin polarization can and has been measured using 
NMR and optical techniques at g = gtr'^^. 

Our calculations are done at zero temperature and in absence of disorder. However, 
disorder is always present in real systems and is crucial to the physics of plateau forma- 
tion in quantum Hall systems. Thermal effects are also expected to be important at the 
temperatures (~ IK) where the experiments are typically done. In fact theoretica l^^i^^ 
and experimental^^ evidence indicates that the skyrmion lattice would have melted at these 
temperatures. It is therefore important to consider the effects of thermal fluctuations and 
disorder before making predictions for experiments. 



A. Spin Polarization at = 

In this section, we present a general proof for the vanishing of the thermal average of the 
spin polarization at (? = in the presence of an arbitrary static disorder potential, Vdis{x, y). 
The classical partition function can be written as, 

gf^E[w,w]+f3 J d?xVdis(x,y)p{x,y) -^^ 

J W,U) 

The spin polarization of the system (of size L) is given by the thermal average 

l^gp^ = A LJl I ~ ^f)E[wM+P S d?^Va,s{x,y)p{x,y) (^52) 

ttL^ 47r Z Jw,w 1 + WW 

Note that at g = 0, the energy functional depends only on the charge density and is 0(3) 
invariant. If we make the transformation of w' 1/w, the charge density and hence 
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the energy density are invariant but the spin density changes sign. Equation ()5.2j) then 
shows that the average spin polarization vanishes. However, this argument for proving that 
spin polarization vanishes on average is technically incorrect since the system is realised in 
a broken symmetry phase. Thus the average in equation (j5.2p is only over configurations 
satisfying the boundary conditions, 

lim w{z,z)=Woo (5.3) 

where i^oo is a fixed value. Since w' and w clearly cannot both satisfy the above condition, 
the above argument is fiawed. 

We now give a more sophisticated argument which takes into account the above point 
and prove the result for a large class of configurations. Without loss of generality, we can 
put Woo = 0. We first consider analytic configurations of an an skyrmion system, 

„=1 [Z - IXn ) 

Where |R"^| < L and the labelling is such that n> m ^ \R^\ > \Rm\- Let R*^ denote the 
zeroes of w. The boundary condition, Woo = 0, ensures that atleast one of the zeroes lies at 
infinity, w can then be written as a product of the positions of the poles and zeros: 



where, 

N 

A = E (5-6) 

n=0 

If A = implies that there is more than one zero at infinity. For every such configuration, 
there are a very large number of nearby configurations of the type A„ = A„ + ^A^, for which 
A 7^ 0. Thus the set of A = configurations has negligible statistical weight. So for a generic 
analytic configuration, only one zero lies at infinity. We concentrate on these configurations 
and rewrite w as, 

A z — R,^^ 

w{z) = - — n , p»«P (5-7) 

Z n=l ^ 

Now consider the configuration. 



X{z- R/) Z-Rn 
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' satisfies the boundary conditions ()5.3p . In the region |-z| << R^, we have 

■pnp N-1 ^ _ -rtnp 



(5.9) 



w{z) 

Thus, for every configuration w{z), there exists another configuration w'{z) which satisfies 
the same boundary conditions and except in a region around one point near the edge, is 
equal to l/w{z). The important fact is that the area of this region does not change as we 
take the thermodynamic hmit, N,L ^ oo, N/{ttL'^) p. For short range interactions we 
then have, 

E[w] = E[w'] + o{L^) 
SP[w] = -SP[w']+oiL'') (5.10) 

Note that a 1/r interaction in two dimensions is short-ranged in this sense. Equation ()5.2|) 
then shows that the average spin polarization density vanishes in the thermodynamic limit. 

We now generalize this proof for arbitrary configurations that have only one zero at 
infinity. Consider configurations that satisfy, 

lim w{z,z) = 

|z|— >oo 

lim zw{z,z) = X (5-11) 

|2:|— >00 



Define, 

w'iz,z)=(^^\ —— (5.12) 



R7_V 1 

z - Rlf I w(z, z) 



It can easily be verified that w' satisfies the boundary conditions ()5.3p and that equations 
(j5.9|) and ()5.10|) hold true. Thus if we assume that, as has been shown for analytic configu- 
rations, the generic configuration has only one zero at infinity, we have proved that the spin 
polarization density will vanish in the thermodynamic limit. 



B. Spin Polarization at gf ~ 

The results obtained in the previous section show that for a system with g exactly zero, 
the spin polarization changes discontinuously from its fully polarized value at z/ = 1 to zero 
when the filling factor changes from 1. In a real system, g will never be exactly equal to 
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zero but will have some very small value. Correspondingly, the single skyrmion at z/ = 1 
will have a large but finite size. The (7 ~ regime discussed in this paper will occur at 
filling factors where this size becomes comparable to the inter-skyrmion spacing, eg. in the 
low g experiments^'^, the single skyrmion spin at z/ = 1 has been estimated to be around 
30-40. As mentioned earlier, for "hedgehog" skyrmions, the spin, s can be related to the 
rms charge radius as < >= 2s (see Appendix 2). Setting the inter-skyrmion spacing, 
a = \J'2,Ti / (1 — u) equal to the charge radius, \/2s yields v ~ 0.9. Thus v < 0.9 will be the 
g ^ Q regime for these skyrmions. 

We have calculated the spin polarisation aX g = O.OS^fQ in both the high and low stiffness 
limits. The skyrmion lattice configuration that minimizes the energy now has a value of A 
slightly away from the symmetric point resulting in a non-zero spin polarisation. The spin 
polarisation is calculated from 



The results are plotted in Figs. 14 and 15. As can be seen they are not very different. The 
spin polarization falls sharply to about 10% of its z/ = 1 value at z/ = .97 and v = .98 in 
the high and low stiffness limits respectively. Both are down to about 1% at z/ = 0.9. At 
this value of g, we have calculated the single skyrmion spin, as described in Appendix 1 to 
be 45 and 33 in the low and high stiffness limits respectively. This implies that the spin 
polarization has fallen to 10% when the charge radius is equal to half the inter-skyrmion 
spacing, i.e. when they are just touching each other. It has fallen to about 1% when the 
radius is roughly equal to the lattice spacing. This is consistent with our general arguments 
of the previous paragraph. 

To summarize, we predict that the spin polarization will fall sharply to about 10% of the 
fully polarized value at z/ = 1 — {71/ 4s), where s is the single skyrmion spin measured at 
u = 1. It will fall to almost zero at around z/ = 1 — (vr/s) and remain zero as z/ decreases 
further. 

These predictions are based on our calculations done at zero temperature and in absence of 
disorder. However previous work indicates that the spin polarization may be fairly insensitive 
to these effects. The Hartree-Fock calculation of Brey et. al.^° was done for the pure system 
at zero temperature. Nevertheless, the calculated spin polarization (at g = go) matched 
fairly well with the experiments of Barret et. al.— which were done at 1.55 K. Thus the 




(5.13) 
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predictions in the previous paragraph may well hold for real systems. 

VI. CONCLUSIONS 

In this paper we have studied the ground state configuration of a lattice of overlapping 
skyrmions in various stiffness regimes. Their behaviour is qualitatively different from small, 
well-separated skyrmions. While well-separated skyrmions can be described well by the 
oriented, charged disc picture, for overlapping large skyrmions this picture breaks down. A 
lattice of large skyrmions is described by the form of the ansatz given in Section III which 
is specified completely in terms of the positions and strengths of the poles and zeros. The 
north poles are the lattice points and the south poles determined by the parameter q. In the 
case of Neel ordered, face- centered rectangular lattices, the south poles lies at the center of 
the unit cell. We then argued that there exists a special value of A = when the charge 
and spin densities are symmetrically distributed around the north and south poles. In the 
case when the lattice was rectangular the charge density formed a smaller face-centered 
rectangular lattice. When the lattice was square the charge density also had the symmetry 
of a smaller square. In all cases the spin density was peaked in opposite directions at the 
north and south poles and the z component of the spin vanishes. Since this distribution 
corresponds to the most symmetric distribution of spin and charge we expect that this also 
minimises the energy. 

We verify the above arguments numerically in Section IV. We find that in the high stiffness 
regime the underlying skyrmion lattice is Neel ordered and rectangular (with aspect ratio 
-\/3) ) corresponding to the charge density having the symmetry of a triangular lattice. The 
value \min which minimises the energy was found to be exactly equal to the value A^ at 
which there was complete symmetry between the north and south poles. In the low stiffness 
regime the underlying lattice and charge density had square symmetry. In the intermediate 
regimes the lattice was rectangular with an aspect ratio between -\/3 and 1. In all cases the 
2;-component of the spin vanished. 

We give a general argument for the vanishing for the spin polarization at = at 
finite temperature and in the presence of disorder. We also discuss some experimental 
consequences. We calculate the spin polarisation at small but non-zero g. We predict that 
this quantity falls sharply to 10% of the fully polarised value at = 1 — 7r/4s and to almost 
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zero aX, V — 1 — n / s. 

We conclude that in the regime of small Zeeman coupling, skyrmions are large and over- 
lapping and the oriented, charged-disc picture breaks down. This is because there is a 
distribution of charge and spin around the south poles also unlike in the charged-disc pic- 
ture where they were distributed only around the north poles. The charge and spin densities 
can have a higher symmetry than the underlying lattice in the overlapping case unlike in 
the well-separated case. The universal features of the classical ground state are (1) Neel 
ordered rectangular lattices of skyrmions with the charge density having the symmetry of a 
face-centered rectangular lattice and (2) vanishing of the 2;-component of the spin. 
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Appendix 1 

In this appendix we calculate the scaling form of the single skyrmion energies with the 
parameter A and also the spin polarisation of a lattice of such skyrmions. The ansatz chosen 
to minimize the energy of a single skyrmion with winding number 1 at = is of the form 

w(z) = -e-^"^ (6.1) 
r 

The ansatz is not analytic anymore and has an exponential damping factor. From the form 
of the individual energies given in the previous sections we can see that they scale with A as 

E,{k,X) = X^E,{k,1) 
Ecouiii^,X) = jEcouiii^A) (6.2) 

From this scaling relation we can calculate the A which minimizes the energy 

3 E,{k, 1) 

Therefore the minimum energy is 

E{k, Xmin) = EgrniK, 1) + XI^^E,{k, 1) + -^E^oul{l^, 1) (6.4) 
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The excess spin per skyrmion corresponding to the minimum configuration is 

2 1 /■ WW 

spin = X^inTT- / T , - „ (6-5) 
zn J 1 + WW 

To extract the spin polarisation of a collection of well separated Nsky number of single 
skyrmions we have to sum over all the individual spins and subtract out the spin of the 

vacuum state {spin^ac)- Thus, 

spiuDL = Ns),y spin - spin^ac 
spinoL spin v 
A ~ Aa 47r 

where A is the area of the sample and A^ is the area of the unit cell of the dilute lattice 
formed by the skyrmions. 



Appendix 2 

In this appendix we show the relation between the spin s of a hedgehog skyrmion and its 
rms charge radius r i.e. (r^) = 2s. Consider a single skyrmion described by the ansatz 

w{z) = ^ (6.7) 
z 

where r is the rms charge radius and F{r) is an analytic function of the rms charge radius. 
The spin of this skyrmion is 

1 r F'^cPr 



Its charge density is 



27rJ r2 + F2 



, . 1 tijdiwdjW 
p[r) — 



2ni (1 + ww) 
J_^j^^wd^ 

2m \l + ww) ^ ■ ' 



Now, 



^'^1 [ wdjW 



(6.8) 



2m J (1 + ww) 
r 2riwd^w 

-TT^ / 7T^ — 6.10) 
2m J {l + ww) ^ ' 
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Using, 



djW 



F'r. 



i.e dj-w — — 



zr 

w 



+ Fdj{l/z) 
& 



dyW 



z 

iw 

z 



(6.11) 



we get 



(6.12) 
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Fig. 1. The face-centered rectangular lattice with basis vectors 

























































— .--.N T 















































































Fig. 2. XY spin configuration in a triangular lattice of well separated skyrmions 

(q= (27r/3,27r/3)). 
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3. Charge density profile in a triangular lattice of well separated skyrmions. The higher 

values are shaded in white. 
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. 4. XY spin configuration in a square lattice of well separated skyrmions (q = (7r,7r)). 
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Fig. 5. Charge density profile in a square lattice of well separated skyrmions. The higher values 

are shaded in white. 





Fig. 6. XY spin configuration in a Neel ordered, rectangular lattice of overlapping skyrmions. 
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Fig. 7. Charge density profile in a Neel ordered, rectangular lattice of overlapping skyrmions. 

The higher values are shaded in white. 




Fig. 8. XY spin configuration in a Neel ordered, square lattice of overlapping skyrmions. 
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Fig. 9. Charge density profile in a Neel ordered, square lattice of overlapping skyrmions. The 

higher values are shaded in white. 



0.004 - 



Fig. 10. Plot of the value of minimzed energy (with respect to A and q). Notice that the 
minimum occurs when the aspect ratio has value \/3. 
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g. 11. Charge density profile in the unit cell of the rectangular lattice (aspect ratio \/3) at 

V = 0.90, in the high stiffness limit. 
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filling factor 



12. Scaling of A with the lattice size paramter e as a function of the filling factor in the high 
stiffness regime. A similar graph is obtained in the low stiffness regime also. 
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Fig. 13. Charge density profile in the square unit ceU in the low stiffness limit at v = 0.98. 




Fig. 14. Plot of the value of the spin polarization as a function of the filling factor in the high 
stiffness regime. The dotted line is the curve calculated from a dilute lattice of isolated skyrmions. 
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Fig. 15. Plot of the value of the spin polarization as a function of the filling faetor in the low 
stiffness regime. The dotted line is the curve calculated from a dilute lattice of isolated skyrmions. 
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